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ABSTRACT 

Cosmic-ray proton and helium spectra have been measured with the balloon-borne 
Cosmic Ray Energetics And Mass experiment flown for 42 days in Antarctica in the 
2004-2005 austral summer season. High-energy cosmic-ray data were collected at an 
average altitude of ~38.5 km with an average atmospheric overburden of ~3.9 g cm -2 . 
Individual elements are clearly separated with a charge resolution of ~0.15 e (in charge 
units) and ~0.2 e for protons and helium nuclei, respectively. The measured spectra at 
the top of the atmosphere are represented by power laws with a spectral index of —2.66 
± 0.02 for protons from 2.5 TeV to 250 TeV and -2.58 ± 0.02 for helium nuclei from 
630 GeV nucleon -1 to 63 TeV nucleon -1 . They are harder than previous measurements 
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at a few tens of GeV nucleoli . The helium flux is higher than that expected from the 
extrapolation of the power law fitted to the lower-energy data. The relative abundance 
of protons to helium nuclei is 9.1 ± 0.5 for the range from 2.5 TeV nucleon" 1 to 63 TeV 
nucleon -1 . This ratio is considerably smaller than the previous measurements at a few 
tens of GeV nucleon -1 . 



Introduction 



Cosmic rays are the product of energetic processes in the universe, and their interactions with 
matter and fields are the source of much of the diffuse gamma-ray, X-ray, and radio emissions 
that are observed. Therefore, the origin of cosmic rays and how they propagate have a major 
impact on our understanding of the universe. Supernova shock waves could provide the power 
required to sustain the galactic cosmic-ray intensity, but details of the acceleration mechanism are 
not completely understood. The shock acceleration mechanism is believed to be a prevalent process 
in astrophysical plasmas on all scales throughout the universe. It has been shown to work in the 
heliosphere, e.g., at planetary bow shocks, at interplanetary shocks in the solar wind, and at the 
solar wind termination shock. 

It is a characteristic of diffusive shock acceleration that the resulting particle energy spectrum 
is much the same for a wide range of shock properties. This energy spectrum, when corrected for 
leakage from the Galaxy, is consistent with the observed spectrum of Galactic cosmic rays. In the 
most commonly used form of t he theory, the characteris tic limiting energy is about Z x 10 14 eV, 
where Z is the particle charge (ILagage &: Cesarskylll983l ). The observed composition should begin 
to change beyond about 10 14 eV, the limiting energy for protons, and the Fe spectrum would start 
to steepen at an energy 26 times higher. In this scenario, protons would be the most dominant 
element at low energies, but heavier elements would become relatively more abundant at higher 
energies, at least up to the acceleration limit for iron. 

Compelling evidence that supernova remnants (SNRs) are common sites for shock acceleration 
of electron s comes from observations of non-thermal synchrotron rad iation from several shell- type 
remnants (jKoyama et al.lll995l : lAllen et al.lll997l : iLeBohec et al.ll2000l ). Non-thermal X-ray spectra 
indicate the presence of very high energy electrons which, at least in the case of SN 1006, have 
energies >2 xlO 14 eV (IKovama et al.lll995l ). These electrons were likely accelerated at the remnant 
because at this energy electrons cannot travel far from their origin before they are attenuated by 
synchrotron losses. T here are other sources of particle acceleration that may also contribute to 
the cosmic-ray beam ( Dermerj boOll ). Recent Chandra X-ray observations of Tycho's SNR have 
shown hot stellar debris keeping pace with an outward-moving shock wave indicated by high-energy 
electrons. 



Semi-direct evidence for t he acceleration of cosmic-ray p rotons could come in the form of 
gamma rays from pion decay (lEllison fc Cassam-Chenail 120051 ). Indeed, the observation of TeV 
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gamma rays, possibly of 7r°-origin, from the SNR RX J1713.7— 3946 ( Enomoto et al.ll2002 ; Aharonian et al. 
20071 ) may have revealed the first specific site where protons are accelerated to energies typical of 
the main cosmic-ray component. Their hadronic origin is yet to be confirmed, but the CANGA- 
ROO collaboration has shown that the energy spectrum of gamma-ray emission from SNR RX 
J1713.7— 3946 matches that expected if the gamma rays are the decay products of neutral pi- 
ons generated in p-p collisions. Although the proton scenario is favored because of the spectral 
shape, gamma rays may originate from either electrons or proto ns. A complete unders tanding of 
gamma-ray emission processes may need a broadband approach (jAharonian et al.ll2006l ). using all 
the available measurements in different wavelength regions. Direct measurements of nuclear parti- 
cle composition changes would provide strong corroborating evidence that shocks associated with 
shell- type SNRs provide the acceleration sites for cosmic rays. 

Shock acceleration is the generally accepted explanation for the characteristic power-law feature 
of cosmic-ray energy spectra, although ground-based measurements have shown that the all-particle 
spectrum extends far beyond the highest energy thought possible for supernova shock acceleration. 
These measurements have also shown that the energy spectrum above 10 16 eV is somewhat steeper 
than the spectrum below 10 eV, which lends credence to the possibility of a different source. Of 
course, the "kn ee" structure might be related to energy-de pendent leakage effects during the propa- 
gation process ( Ptuskin fc Zirakashvililll993l:ISwordv 1995 ) or to other effect s, such as reacceleration 
the galactic wind ( Volk fc Zirakashvililbood ) and acceleration in pulsars ( Bednarek Sz Protheroe 



m 



2002). Whether and how the spectral "knee" is related to the mechanisms of acceleration, propa- 
gation, and confinement are among the major current questions in particle astrophysics. 



2. CREAM Experiment 



The Cosmic Ray Energetics And Mass (CREAM) experiment (jSeo et al.ll2008l ) was designed 
and constructed to extend balloon and space-based direct measurements of cosmic-ray elemental 
spectra to the highest energy possible in a series of balloon flights. The detailed energy dependence 
of elemental spectra at very high energies, where the rigidity-dependent supernova acceleration 
limit could be reflected in composition change, provides a key to understanding the acceleration 
and propagation of cosmic rays. We report in this paper the proton and helium spectra as well as 
their ratios observed from the maiden flight of the CREAM payload in Antarctica. Results from 
the CREAM experimen t such as B/C ratio and heavy elemental spectra are discussed in elsewhere 
(|Ahn et aliboosl . l2009al . hoiol b 



2.1. CREAM Flight 2004 2005 



The first Long Duration Balloon (LDB) flight of the CREAM payload was launched from 
McMurdo Station, Antarctica on December 16, 2004. It subsequently circumnavigated the South 
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Pole three times for a record-breaking duration of 42 days; the flight was terminated on January 
27, 2005. The instrument float altitude remained between 37 and 40 km through most of the flight. 
The corresponding atmospheric overburden was 3.9 ± 0.4 g cm~ 2 . The diurnal altitude variation 
due to the Sun angle change was very small, <1 km, n ear the pole, i.e ., at high latitude, which 
increased as the balloon spiraled out to lower latitudes (|Seo et al.ll2008l ). The temperature of the 
various instrument boxes stayed within the required operational range with daily variation of a few 
°C, consistent with the Sun angle. A total of 60 GB of data including ~^ 
collected. 



10 7 science events were 



The science instrument was supp orted by the command and data module developed by the 
NASA Wallops Flight Facility (WFF) ([Thompson & Stuchliklbood ). This is in contrast to typical 
LDB payloads which utilize the support instrumentation package provided by the Columbia Sci- 
entific Balloon Facility. CREAM was the first LDB mission to transmit all the prime science and 
housekeeping data (up to 85 kbps) in near real-time through the Tracking and Data Relay Satellite 
System (TDRSS) via a high-gain antenna, in addition to having an onboard data archive. To fit the 
data into this bandwidth, science event records excluded information from channels that had levels 
consistent with their pedestal value. This "data sparsification" reduced the average high-energy 
shower event record size by nearly 95%. The science instrument was controlled from a science 
operation center at the University of Maryland throughout the flight after line-of-sight operations 
ended at the launch site. Primary command uplink was via TDRSS, with Iridium serving as backup 
whenever the primary link was unavailable due to schedule or traversing zones of exclusion. The 
nearly continuous availability of command uplink and data downlink allowed rapid response to 
changing conditions on the payload (e.g., altitude-dependent effects) throughout the f light. More 



details about flight ope rations and the data acquisition system are discussed elsewhere (|Yoon et al 
2005al : lzinn et al.ll2005h . 



2.2. CREAM Instrument 



The instrument was designed to meet the challenging and conflicting requirements to have a 
large enough geometry factor to collect adequate statistics for the low flux of high-energy particles, 
and yet stay within the weight limit for near-space flights (jAhn et al.ll2007al ). It was comprised of 
a suite of particle detectors to determine the charge and energy of the very high energy particles. 
As shown schematically in Figure (U the detector configuration included a timing charge detector 
(TCD), a transition radiation detector (TRD) with a Cerenkov detector (CD), a silicon charge 
detector (SCD), hodoscopes (HDS), and a tungsten/scintillating fiber calorimeter. Starting from 
the top, the TCP c onsists of two crossed layers of four 5 mm thick and 1.2 m long plastic scintillators 
(jAhn et al.ll2009bl ). It defines the 2.2 m 2 sr trigger geometry and determines charge based on the 
fact that the incident particle enters the TCD before developing a shower in the calorimeter, and 
the backscattered albedo particles arrive several nanoseconds later. A layer of scintillating fibers, 
S3, located between the carbon target and the tungsten calorimeter provides a reference time. 
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Fig. 1. — CREAM detector configuration. From the top are shown the TCD, upper TRD1, CD, 
lower TRD2, SCD, Hodoscope SO/SI, carbon targets, S2, S3 and the calorimeter. 



The TRD determines the Lorentz factor for Z ^ 3 nuclei by measuring transition X-rays 
using thin- wall gas tubes. Transition radiation is produced when a relativistic particle traverses 
an inhomogeneous medium, in particular the boundary between materials of different dielectric 
properties. The TRD consists of a foam radiator and 16 layers o f proportional tubes filled with a 
mixture of xenon (95%) and methane (5%) gas (jAhn et al.l l2008). The CD between the two TRD 
sections provides low-energy particle rejection at the flight site, Antarctica, where the geomagnetic 
cutoff is low. It also provides additional charge identification. 



The SCD is comprised of 380 fim thick Si sensors (|Park et al.ll2007al ). It is segmented into 
pixels, each about 2.12 cm 2 in area to minimize multiple hits in a segment due to backscattered 
particles. The targets are comprised of blocks of densified graphite cemented in carbon/epoxy 
composite cradles. The vertical thickness of the carbon targets is about 0.5 interaction lengths. 
They force hadronic interactions in the calorimeter, which measures the shower energy and provides 
tracking infor mation to deter mine which segment (s) of the charge detectors to use for the charge 
measurement (jSeo et al.lll996l ). 



The calorimeter consists of 20 tungsten layers interleaved with scintillating fiber ribbon layers 
which are alternately oriented in the x- and y-directions. Each tungsten layer is 1 radiation length 
thick to sample the shower every radiation length. Each layer consists of fifty 1 cm wide and 
0.5 mm thick fiber ribbons to measure the longitudinal and lateral distributions of the shower. 
The light signal from each ribbon is collected by means of an acrylic light-mixer coupled to a 
bundle of clear fibers. This is split into three sub-bundles, each feeding a pixel of a hybrid photo 
diode (HPD). In this way the wide dynamic range of the calorimeter is divided into three sub- 
ranges (low, mid, high) with different gains, chosen to match the dynamic range of the front-end 
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electronics (|Lee et al.l 120061 ) . Tracking for showers is accomplished by extrapolating each shower 
axis back to the charge detectors. The HDS SO/SI and S2, comprised of 2 mm thick and 2 mm wide 



scintil l ating fibers, provide ad ditional tracking information above the tungsten stack (|Yoon et al 



( Ahnetal 



2001 



2005b; Marrocch esi et al.ll2004l ). The tracking uncertainty is smaller than the pixel size of the SCD 



Tracking for non- interacting particles is achieved in the TRD with better accuracy (1 mm res- 
olution with 67 cm lever arm, 0.0015 radians). The TRD and calorimeter have different systematic 
biases in determining particle energy. The use of both instruments allows in-flight cross-calibration 
of the two techniques and, consequently, provides a powerful method for measuring cosmic-ray en- 



(Lee et al. 


2006 


: Park et al. 


2004; 


Ahn et al. 


2007a) 



3. Data Analysis 

The main trigger conditions for science events were (1) significant energy deposit in the 
calorimeter for high-energy particles or (2) large pulse height, Z > 2, in the TCD for heavy nuclei. 
The former requires each of six consecutive layers in the calorimeter to have at least one ribbon 
recording a deposit of more than 45 MeV. The high-energy shower events that meet this calorimeter 
trigger condition were used in this analysis. 



3.1. Event Selection 

The ribbon with the highest energy deposit and the neighboring ribbons on both sides were 
used to determine the position in each layer of maximum energy deposits. The shower axis was 
reconstructed by a least-sq uares fit of a strai ght line through a combination of these hit positions 



in the XZ and YZ planes (lAhn et al.ll2007bl ). Hits not along the straight line were excluded from 
the fit. The resulting trajectory resolution is ~1 cm when projected to the SCD. The reconstructed 
trajectories were required to traverse the SCD active area and the bottom of the calorimeter active 
area. 

At this stage non-interacting particles are removed, but some events have their first hadronic 
interaction in the calorimeter layers instead of the carbon targets. These late interacting events 
could result in an underestimation of deposited energy, or misidentification of charge due to large 
uncertainties in the trajectory reconstruction. Since their longitudinal shower profiles are different, 
events with small energy deposit in the top few layers of the calorimeter were removed to ensure 
that the selected events had their first interactions either in the carbon targets or in the top of the 
calorimeter. 
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3.2. Charge Determination 

In order to determine the incident particle charge, the reconstructed shower axis from the 
calorimeter was extrapolated to the SCD and a 7 x 7 pixel area, about 10 x 10 cm 2 , centered on 
the extrapolated position, was scanned to seek for the highest pixel signal. The scanned area was 
optimized to sustain the charge identification efficiency of 99% in all energy bins, accounting for 
dead and noisy SCD and calorimeter channels (~15% and 13%, respectively), and determined to 
be a 7 x 7 pixel area. That highest pixel signal was then corrected for the particle path length 
(calculated from the reconstructed incidence angle) in the sensor. The signal reflects the ionization 
energy loss per unit path length {dE/dx) of an incident particle in the SCD. The energy loss is 
proportional to Z 2 . According to Monte Carlo (MC) simulations and beam tests, the expected 
contamination from secondary particles back scattering f rom the calorimeter is <3% when this 



tracking-based selection method is used (jPark et al.ll2007bl ). The resulting SCD signal distribution 



is shown in Figure [2j Events with Z < 1.7 were selected as protons, while events with 1.7 ^ Z 
< 2.7 were selected as helium nuclei. The charge resolutions are estimated as ~0.15 e and ~0.2 e 
for protons and helium nuclei, respectively. The proton and helium losses due to dE/dx Landau 
tails were corrected by charge selection efficiencies, which will be discussed in Section 13.61 The 
proton events in the helium range were removed as a background in the helium selection, and the 
helium events in the proton range were removed as a background in the proton selection, which 
will be discussed in Section T3.5I Unstable SCD channels identified by their large root-mean-square 
pedestal variations throughout the flight were excluded from the analysis. Including dead or noisy 
channels, ~15% of the total 2,912 SCD channels were masked. 



3.3. Energy Measurement 

An ionization calorimeter is the only practical way to measure the energy of protons and 
helium nuclei above ~1 TeV, but calorimeters with full containment of ha dronic showers are too 



massive to be incorporated into space-based or balloon-borne experiments (jGanel et al.lll999l ). A 
thin calorimeter offers a practical approach but the calorimeter calibration requires the use of 
accelerator beam particles having known energy. The CREAM calorimeter was calibrated before 
the flight with electron beams at the European Organization for Nuclear Research (CERN). Each 
of the 1000 fiber ribbons was exposed to 150 GeV electrons. The responses from the 50 ribbons in 
a given layer are equalized by moving the detector in steps of 1 cm vertically or 1 cm horizontally, 
so the electron beam is centered each time on the center of a different ribbon in each X or Y layer. 

The calorimeter was designed to measure the energy deposit from showers initiated by nuclei 
with energies up to 10 15 eV and higher. Its sampling fraction for isotropically incident TeV proton 
showers initiated in the graphite targets is about 0.13% of the parent's energy in the active media. 
With electron test beam energies of 150 GeV or less, only 8—10 layers around the shower maximum 
register enough scintillation to allow calibration. To address this, the calibration scan was carried 
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Fig. 2. — Charge distribution of selected events determined with the SCD. The number of protons 
and helium nuclei are identified from the charge distribution in each energy bin. 



out in three sets of runs by ex posing the c alorim eter from the bottom with additional targets along 
the beam line as described in 

The energy deposit expected along the shower core in each layer was calculated using MC 
simulations of electron showers. Conversion factors from analog-to-digital conversion (ADC) signals 
to MeV were obtained from the ratio of MC simulation of the energy deposited in each ribbon to 
the measured ADC sig nal from the cali b ration beam test . The MC simulations were based on 
GEANT/FLUKA 3.21 ferun et al.lll984l : iFasso et al.lll993l '). The ADC signals were corrected for 
the HPD quantum efficiency and gain difference from the different HPD high-voltage settings 
between the beam test and the flight. 

Inter-calibration between the low- and mid- energy ranges, and between the mid- and high- 
energy ranges were carried out with flight data by comparing the signals from two ranges of the same 
ribbon generated by the same shower. None of the proton and helium event candidates saturated in 



Yoon et al. 


(2005c. 


2007 


) and 


Ahn et al. 



The calorimeter, HDS, and SCD were also exposed to n uclear fragments ( A/Z = 2) of a 158 



GeV nucleon 1 Indium beam at CERN ( Yoon et al. 



2004 



2007 



Ahn et al. 



2006a: 



Marrocchesi et al 



Park et al.ll2004l ). The energy response was linear up to the maximum beam energy of ~9 
TeV. Above the available accelerator beam energy, MC simulations indicate that the calorimeter 
response is quite linear in the CREAM measurement energy range. Simulations also indicate that 
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the calorimeter energy resolution is nearly energy independent ( Ahn et aOhoOll ). Nevertheless, our 
energy deconvolution i ncluded correction s for the small energy dependence of the energy resolution 
due to shower leakage (|Ahn et al.ll2009al ). 



3.4. Spectral Deconvolution 



Entries in the deposited energy bins were deconvolved into incident energy bins using matrix 
dations. The counts, Ni nCt i, in incident energy bin i were estimated from the measu red counts, 
dep,jj hi deposited energy bin j by the relation (jBucklev et al. 1994 ; Ahn et al. 2006b ) 



,3i 



(1) 



where matrix element Pu is a probability that the events in the deposited energy bin j are from 
incident energy bin i. The matrix element Py was estimated from the response matrix generated 
by MC simulation results obtained separately for protons and helium nuclei. The response matrix 
and corresponding deconvolution matrix were generated and tested by varying the indices between 
—2.5 and —2.8. We verified that the flux deconvolution process was not sensitive to the assumed 
spectral index used, within that range, to generate the matrix elements. 



The MC simu l ations for helium and heavy nuclei used FRITIOF /RQMD /DPM JET-II (|Kim et al 



1999 



Wang et al 



2001 



int erface d to the GEANT/FLUKA 3.21 hadronic simulation package. 



FRITIOF (jAndersson et al.1 Il993l ) is based upon semiclassical considerations of string dynamics 
for high-energy hadronic collisions. The relativistic quantum molecular dynamics (RQMD) model 
was adopted for simulations of heavy ions for energies in the center-of-mass frame less than 5 
GeV nucleon -1 . RQMD is a semi classical mic roscopic approa ch which combines classical prop aga- 



tion with stochastic interactions (|Sorgdll995l ). DPM JET-II (jRanft 



1995 



Ferrari et al.lll996l ) was 



based on the dual parton model, a framework for hadron-hadron interactions and production in 
hadron-nucleus and nucleus-nucleus collisions at high energies. 



3.5. Background Corrections 

The primary background is comprised of events with misidentified charge, which result mainly 
from secondary particles generated by interactions above the SCD or from particles back-scattered 
from the calorimeter. This is the case for the protons; however, there is an additional cause of 
misidentified events for helium nuclei: the proton dE/dx Landau tail. Misidentified event counts of 
protons and helium nuclei were estimated from the MC simulations with a power-law input spec- 
trum. Due to the Landau tails, back-scattered and secondary particles, 5.1% of measured protons 
were misidentified helium nuclei and 6.8% of measured helium nuclei were misidentified protons, 
as shown in Table [TJ About 0.2% of incident carbon nuclei were identified as protons, and 2.8% of 
incident carbon nuclei were misidentified as helium nuclei, using the energy spectra of individual 
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cosmic-rays compiled by IWieb el- Sooth et al.l (|1998l ). Less than 1% of trigger and reconstructed 
protons and helium events are from secondary particles. Additional background comes from the 
events that are not within the geometry, but which satisfy the trigger and reconstruction conditions; 
they are either entering the instrument acceptance from outside the SCD area or exiting the side 
of the calorimeter instead of the bottom. According to MC simulations, this is about 3.6% and 
4.0% of the selected events for protons and helium nuclei, respectively. The total background was 
9% for protons and 11% for helium nuclei. 



3.6. Absolute Flux 

The measured spectra are corrected for the instrument acceptance as shown below to obtain 
the absolute flux F: 

F = ^ I (2) 

dE GF e T rj' v ; 

where dN is the number of events in an energy bin, dE is the energy bin size, GF is the geometry 
factor, e is the efficiency (defined below), T is the live time, and r] is the survival fraction after 
accounting for atmospheric attenuation. The geometry factor was calculated to be 0.43 m 2 sr using 
an MC simulations by requiring the extrapolated calorimeter trajectory of the incident particle to 
traverse the SCD active area and the bottom of the calorimeter. Out of 42 days of the flight, the 
stable period was about 24 days when no commands were sent, e.g., for instrument tuning, power- 
cycle, or high-voltage adjustments. After the dead-time correction, the live time, T, of 1,099,760 s 
was used for this analysis. 

Efficiency. The efficiency, e in Equation ([2]) includes efficiencies from all analysis steps, includ- 
ing trigger condition, event reconstruction, charge identification, and removing events with late 
interactions: 

£ — E-trig &rec &sel ^-charge • (3) 

The trigger efficiency, £t r i g , was obtained from the fraction of events satisfying the trigger condition 
among all events within the geometry, i.e., passing through the bottom of the calorimeter and the 
SCD active area, using MC simulations. This is energy dependent at low energies where the trigger 
is not fully efficient. Above 3 TeV, it is nearly constant around 76% for protons and 91% for 
helium nuclei, respectively. The reconstruction efficiency, e rec , was taken to be the ratio of events 
satisfying the reconstruction and trigger conditions to events satisfying only the trigger condition. 
The reconstruction efficiency was 98% for protons and 99% for helium nuclei, respectively, based 
on MC simulations. The event selection efficiency, e se i, was estimated with the MC simulations 
after removing events with late interactions and was 90% protons and 96% for helium nuclei. The 
charge efficiency, £ c harge, takes into account lost events due to the noisy or dead SCD channels, 
interactions above SCD and misidentified charges. It was calculated to be 77% for protons and 67% 
for helium nuclei, respectively, using MC simulations. The efficiencies are summarized in Table [TJ 



The trigger efficiency for proton and helium nuclei cannot be estimated with the flight data, 
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since we do not know how many un-triggered events occurred. However, the event selection ef- 
ficiency, s se i, and charge efficiency, s charge, were estimated in a limited way using flight data for 
combined protons and helium nuclei events. It is not as accurate as individual MC simulations 
because the composition (abundance) of the incident particles is unknown. When the abundance 
ratio of protons and helium nuclei was assumed to be 1:1 and the abundance of heavy nuclei above 
helium nuclei was ignored, the combined efficiencies were 68% from the flight data and 67% for the 
MC simulations. 

Interactions in air. The attenuation loss due to the atmospheric overburden, 3.9 ± 0.4 g cm -2 , 
was corrected for survival fractions of protons and helium nuclei. This air depth was measured 
by pressure sensors during the flight. Interaction cross sections have been measured in many 



fixed target experiments, and cross sec tions are known up to a few tens of GeV (IHaeen et a 



1977 



(1977) 



Webber et al.lll990l ; iPapini et al.lll996l ). We used the cross section formula from lHagen et al. 
to calculate interaction lengths and survival fractions for protons and helium nuclei. The mean 
incident angle of 35°, estimated from the flight data, was used to estimate the losses. The survival 
fraction, n, used to characterize atmospheric attenuation was determined to be 95% for protons 
and 91% for helium nuclei, respectively. 



The ratio of secondar y to primary protons and helium nucle i in the atmosphere above GeV 



energies has been reported (IKawamura et al 



1989 



Abe et al.ll2003h . IPapini et all (j 1996h calculated 



that the secondary to primary proton ratio at an air depth of 3 g cm -2 was less than 1% above 
40 GeV, and the secondary to primary helium nuclei ratio was less than 2% at 10 GeV nucleon -1 . 
Our MC simulations showed that the fraction of secondary protons and helium nuclei produced 
from carbon and iron nuclei interactions in the air was less than 1% at 10 TeV. 

Energy-bin representation. For the number of events (dN) in each energy bin with upper- and 
lower-energy limits, Ej+i and Ej, respectively (dE = Ej + \ — Ej), the differential flux is dN/dE at 
E m , where E m can be taken as the arithmetic mean of Ej and Ej + \ in logarithmic range or else 
using a suitably weighted average of Ej and -Ej+i • We als o investigated an alternative procedure 
to determine E m , as suggested by lLafferty fe Wyattl (|1995l ): 

'■Ej4 



1 



E 



'j+i 



■'3+1 



f(E)dE. 



(4) 



Table 1. Efficiencies and backgrounds for absolute flux 



Efficiency and Background 


Proton (%) 


Helium (%) 


Trigger efficiency 


76 ± 2 


91 ± 1 


Reconstruction efficiency 


98 ± 1 


99 ± 1 


Late interaction events efficiency 


90 ± 1 


96 ± 1 


Charge selection efficiency 


77 ± 2 


67 ± 2 


Background from reconstruction 


3.6 ± 0.1 


4.0 ± 0.2 


Background from misidentified charge 


5.1 ± 0.2 


6.8 ± 0.2 
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For a power-law spectrum, f(E) = AE 7 , E m can be calculated as, 



E*77 - E)-^ v -i/7 

^=( , P 3+1 ^ • (5) 



.(^• +1 -^)(l-7)' 

In this analysis, i£ m was used and the difference between E m and the center of the bin in logarithmic 
range is less than 1%. 



3.7. Uncertainties 



The statistical uncertainty in each energy bin was estimated b y the relation 8Nj nCi i = 5{Y^ - Pij^depj): 



considering 68.3% the Poisson confidence interval determined bv lFeldman <fe Cousins! (|1998l ). The 
uncertainties were estimated by propagating uncertainties of measured entries in each bin and un- 
certainties of dec onvolution components, Pij, from MC simulations, while in the paper reported by 



Ahn et al.l (|2010l ) , uncertainties were estimated by propagating uncertainties from measured entries 
with P^. This estimation gives more conservative results than the reported results. 

Several sources of systematic uncertainties were identified. The systematic uncertainties for 
efficiencies and backgrounds were estimated within each energy range to account for the energy- 
dependent effects determined using MC simulations. They are summarized in Table [TJ Efficiency 
uncertainties were about 1-2% and background uncertainties were about 5%. 

The geometry factor uncertainty was 2% for both protons and helium nuclei; it was esti- 
mated with MC simulations. The precision of estimated live-time fraction was about 3.3% and 
the accuracy of estimated dead time due to timeouts in TCD readout, which delayed processing, 
was about 2.6%. The overall uncertainties for the estimated live time were 4% for both protons 
and helium nuclei. The systematic uncertainties for the survival fractions in the atmosphere were 
calculated analytically. The p-p cross section difference between 10 TeV a nd 100 TeV is abou t 



28%, according to the most recent reference from the Particle Data Group (jAmsler et al.ll2008l ). 
Using a conservative estimate of 30% for cross section uncertainties, the estimated uncertainties 
of survival fractions were 2% and 3% for protons and helium nuclei, respectively. The range of 
incident angle was from 0° to 66°. The uncertainty in correcting for atmospheric losses introduced 
by using an assumed mean incident angle was at the level of 1% for protons and 1.6% for helium 
nuclei. The energy calibration accuracy was found to be 1%. The systematic uncertainties of the 
measured number in each energy bin, considering the 1% energy calibration accuracy, were 3% 
for both protons and helium nuclei. To estimate uncertainties in the spectral deconvolution, the 
unfolding procedure was repeated by varying input spectral indices. The difference of the proton 
fluxes varying the input spectral indices between 2.64 and 2.68 was less than 1%. Similarly, the 
helium flux difference was also less than 1% for input spectra between indices 2.56 and 2.60. The 
overall systematic uncertainties were found to be 9% for both protons and helium nuclei. These 
systematic uncertainties are energy independent. They do not change the spectral shape, but they 
might shift the normalization of the spectra up or down. 
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Results 



The measured proton fluxes from 2.5 TeV to 250 TeV and helium fluxes from 630 GeV 
nucleon -1 to 63 TeV nucleon" 1 at the to p of the atmosph ere are given in Tables [2] and [3l while 
previously reported results in the paper bv lAhn et al.l (|2010l ) are presented in a plot. The statistical 
uncertainties were re-estimated, as discussed in Section T3.71 The CREAM proton and helium spec- 
tra are each consistent with a single power law over the measured range. The best-fit parameters 
for the spectra for protons and helium nuclei are represented by 



^5- = <fr E 13 (m 2 sr s GeV nucleon 1 ) 1 . 
dE 



(6) 



The best-fit parameters for the spectra for protons and helium nuclei are given by (|Ahn et al.ll2010l ): 

(7.8 ± 1.9) x 10 3 (m 2 srs)~ 1 (GeV nucleon" 1 ) 1 - 66 , 
P P = 2.66 ± 0.02, 



0-P 



(7) 
(8) 



and 



$ He = (4.2 ±0.8) x 10 2 (n^srs)" 1 (GeVnucleon" 1 ) 
(3 He = 2.58 ± 0.02. 



1a 1.58 



(9) 
(10) 



The spectral indices for proton and helium nuclei were calculated both with the least squares 
fit and maximum likelihood method. The results from both methods were consistent. Uncertainties 
for the spectral indices were estimated with the maximum likelihood method. 

The CREAM proton spe ctrum is harder than previous measureme nts at lower energi es such 
as AMS dAguilar et alJ I2OO2I ). 2.78 ± 0.009 at 10-200 GV and BESS (jHaino et alJboO^h . 2.732 
± 0.011 from 30 GeV to a few hundred GeV. Likewise, the CREAM helium spectrum is harder 
than AMS, 2.740 ± 0.01 at 20-200 GV and BESS, 2.699 ± 0.040 from 20 GeV nucleon" 1 to a few 
hundred GeV nucleon" 1 . 

Table 2. Proton differential flux measured with CREAM 



Energy Bin Ran 


ge 


Flux ± Stat. 


(GeV) 




(m 2 sr s GeV)" 1 


2.5 x 10 3 


- 4.0 x 10 3 


(3.72 ±0.10) x 10- 6 


4.0 x 10 3 


- 6.3 x 


10 3 


(1.10 ±0.04) x 10~ 6 


6.3 x 10 3 


- 1.0 x 


10 4 


(3.19 ±0.19) x 10" 7 


1.0 x 10 4 


- 1.6 x 


10 4 


(9.47 ±0.80) x 10" 8 


1.6 x 10 4 


- 2.5 x 


10 4 


(2.80 ±0.35) x 10~ 8 


2.5 x 10 4 


- 4.0 x 


10 4 


(8.1 ± 1.5) x 10~ 9 


4.0 x 10 4 


- 6.3 x 10 4 


(2.2 ±0.6) x 10" 9 


6.3 x 10 4 


- 1.0 x 10 5 


(6.ll 2 ; 6 2 ) x 10-w 


1.0 x 10 5 


- 1.6 x 


10 5 


(1.8lJ'2) x 10- 10 


1.6 x 10 5 


- 2.5 x 10 5 


(4.2l^ 4 ) x 10- 11 
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Figure [3] compa res our me asured spectra with previous m easurements: AMS, BE SS, CAPRICE98 



Boezio et allbood) ATIC-2 (jPanov et al.1 bood l . JACEE (|Asakimori et alJll998h and RUN JOB 



Derbina et al.ll2005l ). The error bars shown in the figures represent the statistical uncertainties. 



The CREAM results are consistent with JACEE where its measurement energy range overlaps with 
CREAM but indicate higher fluxes, particularly for helium, with respect to RUNJOB. The proton 
and helium fluxes are both higher than that expected by extrapolating the power law fitted to the 
lower-energy measurements, which verifies that our TeV spectra are harder than the lower-energy 
spectra. At 20 TeV nucleon -1 the helium flux measured by CREAM is about 4cr higher than the 
flux expected from a power-law extrapolation of the AMS helium flux and spectral index. 

The proton to helium ratio as a function of energy provides insight into whether the proton 
and helium spectra have the same spectral index. This has long been a tantalizing question, mainly 
because of the limited energy range individual experiments could cover. The ratio from the first 
CREAM flight provides a much needed higher energy, low-statistical uncertainty, measurement 



The r atio is compared with pr evious measurements in Figure [H ATIC-2, CAPR ICE94 (jBoezio et al 



19991 ). CAPRICE98, JACEE (jAsakimori et alJll993U ). LEAP (|Seo et al.lll99ll ). and RUNJOB. The 
CREAM ratios are consistent with JACEE where its measurement energy range overlaps. The 
measured CREAM ratio at the top of the atmosphere is on average 9.1 db 0.5 for the range from 
2.5 TeV nucleon -1 to 63 TeV nucleon" 1 , which is significantly lower than the ratio of ~20 obtained 
from the lower-energy measurements. 



5. Discussion and Conclusion 



The energy spectra of primary cosmic rays are known with good precision up to energies around 
10 11 eV, where magnetic spectrometers have been able to carry out such measurements. Above 
this energy the composition and ene rgy spectra are not accurately known, a l though there have 
been some pioneering measurements (jMiiller et al.lll99ll : lAsakimori et al.lll998l : lApanasenko et al 



Table 3. Helium differential flux measured with CREAM 



Energy Bin Range 
(GeV nucleon -1 ) 



Flux ± Stat. 
(m 2 sr s GeV nucleon -1 ) 



6.3 x 10 2 


- 1.0 x 10 3 


(1.42 ±0.04) x 10~ 5 


1.0 x 10 3 


- 1.6 x 10 3 


(4.35 ±0.16) x 10 -e 


1.6 x 10 3 


- 2.5 x 10 3 


(1.31 ±0.07) x 10 -6 


2.5 x 10 3 


- 4.0 x 10 3 


(3.83 ±0.31) x 10 -7 


4.0 x 10 3 


- 6.3 x 10 3 


(1.27 ±0.14) x 10 -7 


6.3 x 10 3 


- 1.0 x 10 4 


(4.19 ±0.64) x 10 -8 


1.0 x 10 4 


- 1.6 x 10 4 


(1.15 ±0.27) x 10 -8 


1.6 x 10 4 


- 2.5 x 10 4 


(3.4±i;J) x 10-s 


2.5 x 10 4 


- 4.0 x 10 4 


(8.2+ 4 ^) x 10 -10 


4.0 x 10 4 


- 6.3 x 10 4 
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Fig. 3. — CREAM Proton and helium differential Flux ■ E 2 ' 75 in GeV nucleon -1 at the top of 
the atmosphere. The CREAM proton and helium spectra (filled circles) are shown together with 
previous measurements: BESS (squares), CAPRICE98 (downward triangles), AMS (open circles), 
ATIC-2 (diamonds), JACEE (stars) and RUN JOB (crosses). The lines represent power-law fits 
with spectral indices of -2.66 ± 0.02 for protons and -2.58 ± 0.02 for helium nuclei, respectively. 



2001). The collecting power of CREAM is about a factor of two larger than that of ATIC for 
protons and helium nuclei and, considering the much larger geometry factor of the TRD, about a 
factor of 10 larger for heavier nuclei. TRACER has a larger geometry factor than CREAM, but a 
smaller dynamic charge range {Z = 8-26) was reported for its 10-day Antarctic flight. Although 
its dynamic charge range was improved to Z =3-26 for its ~4 day flight from Sweden to Canada 
in 2006, it is still insensitive to protons and helium nuclei. 

The CREAM payload maintained a high altitude, corresponding to an atmospheric overburden 
of 3.9 g cm -2 for vertically incident particles. That implies about 6.8 g cm -2 at the maximum 
acceptance angle for this analysis, which is smallest among comparable experiments. For example, 
the average vertical depth for RUNJOB was more than twice that of CREAM, due to its low 
flight altitude. Considering the RUNJOB acceptance of particles at large zenith angles, its effective 
atmospheric depth was as large as 50 g cm -2 . For that depth, large corrections are required to 
account for the fact that 41% of protons would have interacted before reaching the detector. 
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Fig. 4. — Ratio of protons to helium nuclei as a function of energy in GeV nucleon . The 
CREAM (filled circles) ratio of proton to helium is compared with previous measurements: ATIC- 
2 (diamonds), CAPRICE94 (upward triangles), CAPRICE98 (downward triangles), LEAP (open 
circles), JACEE (stars), and RUN JOB (crosses). 



The CREAM calorimeter is much deeper than either that of JACEE or RUNJOB, so it pro- 
vides better energy measurements. CREAM also has excellent charge resolution, sufficient to clearly 
identify individual nuclei, whereas JACEE and RUNJOB reported elementa l groups. Our o b serva - 
tion did not confirm a soft er spectrum of protons above 2 TeV reported by iGrigorov et al.l (|1970l ) 



or a bend around 40 TeV (jAsakimori et al.lll993al ). An increase in the flux of helium relative to 



protons could be inte rpreted as eviden ce for two different types of sources for protons and helium 
nuclei as proposed by iBiermannl (j 19931 ). The observed harder spectra compared to prior low-energy 
measurements may require a significant modification of conventional acceleration and propagation 
models, with significant impact for the interpretation of other experimental observations. 

The CREAM experiment was planned for Ultra Long Duration Balloon (ULDB) flights lasting 
about 100 days with super-pressure balloons. While waiting for development of these exceptionally 
long flights, the CREAM instrument has flown five times on LDB flights in Antarctica. It should 
be noted that a 7 million cubic foot (~0.2 million cubic meters) super-pressure balloon was flown 
successfully for 54 days during the 2008-2009 austral summer season. As ULDB flights become 
available for large science payloads, long-duration exposures can be achieved faster and more effi- 
ciently with reduced payload refurbishment and launch efforts. Whatever the flight duration, data 
from each flight reduces the statistical uncertainties and extends the reach of measurements to 
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energies higher than previously possible. 
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